Dephasing of Atomic Tunneling by Nuclear Quadrupoles 
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Recent experiments revealed a most surprising magnetic-field dependence of coherent echoes in 
amorphous solids. We show that a novel dephasing mechanism involving nuclear quadrupole mo- 
ments is the origin of the observed magnetic-field dependence. 
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Until recently it was the general believe that the di- 
electric properties of insulating glasses - free of magnetic 
impurities - are largely independent of external magnetic 
fields. New investigations, however, have shown that the 
dielectric properties of certain multi-component glasses 
at very low temperatures are strongly influence by a mag- 
netic field I |, I, |, |, §. In particular, the low- 
frequency dielectric susceptibility and the amplitude of 
spontaneous polarization echoes generated in these amor- 
phous materials show a striking non-monotonic depen- 
dence on applied magnetic field. 

Since the low-temperature properties of glasses are 
governed by atomic tunneling systems (for recent reviews 
see [|[ it has been speculated whether and how a mag- 
netic field can couple to quantum tunneling. Two models 
have been proposed, that relate the magnetic-field depen- 
dence to the Ahanorov-Bohm phase of a charged particle 
moving along a closed loop Very recent polar- 

ization echo experiments, however, indicate that such a 
periodic variation of the tunnel splitting is not the origin 
of the observed magnetic field effects [|[ 0. In contrast, 
these experiments strongly suggest that nuclear magnetic 
moments play a crucial in the observed anomalies. 

In this paper we discuss how nuclear magnetic and 
quadrupolar moments influence atomic tunnel states. Af- 
ter a brief reminder of two-pulse echoes of two-level sys- 
tems and the nuclear spin hamiltonian, we give the echo- 
amplitude correction factor due to nuclear spins and eval- 
uate this expression for the limiting cases of weak and 
strong magnetic fields. Finally, we compare our theory 
with recent data for several glasses. 

Two-level systems (TLS) in glasses arise from double- 
well potentials with asymmetry A and tunnel matrix el- 
ement Aq, 



(1) 



where cr^ = ±1 is the reduced two-state variable that in- 
dicates the states localized in the two wells ^, ^. The 
eigenstates tpo and tpi are separated by the energy split- 
ting E = -\/A§ + A2 (see Fig. 1.) In a two-pulse echo 
experiment, the first pulse creates a coherent superpo- 
sition of the ground state and excited state, with a rel- 
ative phase factor ^-^Et/h^ Because of the dispersion of 




FIG. 1: Two-level system with energy splitting E, eigenfunc- 
tions -f/jo and ■i/'i , and the corresponding quadrupole quantiza- 
tion axes uo and ui. 



the two-level splitting E, the corresponding macroscopic 
polarization decays rapidly. After a waiting time t^, 
the second pulse exchanges the amplitudes of these two 
states; the resulting phase factor Q-^E(t-t^)/h jgg^fjg 
revival of the coherent polarization, and the "echo" is 
observed at a time t after the second pulse. 



Po{t, t^) = ^ cos [uji{t - t^)] , 



(2) 



where the sum runs over all TLS with tunnel frequency 
LUi — Ei/h and effective dipole moment Ai. 

Such a tunnel system involves several atoms, each of 
which may carry a nuclear magnetic dipole and an elec- 
tric quadrupole. For the sake of simplicity, we consider a 
single atom whose nucleus is in a state of total angular 
momentum I where = h^I{I + 1). This "nuclear spin" 
results in a magnetic moment gfi^l/h, with the Lande 
factor g and the nuclear magneton /in = 5 x 10^^^ J/T. 
In the case / > 1, the orbital motion of the protons is 
related to an electric quadrupole moment p^ ; for a nu- 
clear charge distribution p{r) oriented along the axis e 
one finds 



Q = 



[3( 



r • e 



P(r). 



(3) 
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The magnetic dipole couples to the external field B = 
Bbz and the quadrupole moment to the electric field gra- 
dient (EFG) that is given by the curvature of the crystal 
field potential <j>{r). We consider the simplest case of a 
single diagonal term (f>" = (u • V)^0 along the axis u. 
Then the spin hamiltonian reads as p2[ 



and eigenfunctions in the upper and lower tunnel states, 



h"Q ill - /(/ + 1) 

4 /(2/-1) 



(4) 



with the projections of the nuclear spin operator on the 
axes defined by the EFG, /„ = (u • l)/fi, and the mag- 
netic field, Jz — Iz/fi- 

For zero magnetic field, u is the appropriate quanti- 
zation axis; with — rr? and m = the hamil- 
tonian is diagonal. (/ = 1 and 1 = 3/2 give rise to a 
doublet, 1 = 2 and / = 5/2 to a triplet, etc.) In the 
opposite case of zero EFG, the usual choice 1^ = m gives 
the (2/ -I- 1) Zeeman levels mgii-^^B. In general, the axes 
defined by the magnetic field, e^, and the EFG, u, are 
not parallel, i.e., the operators /„ and 1^ cannot be diag- 
onalized simultaneously, thus resulting in a more compli- 
cated situation if both the magnetic field and the EFG 
are finite. 

For asymmetric TLS, one of the eigenstates, say the 
ground state Voj has a large probability amplitude in 
the left well, whereas the excited one, Vii has a larger 
amplitude in the right well (see Fig. 1.) In an amor- 
phous or disordered solid, the crystal field, and thus the 
quadrupole quantization axis, are not the same in the 
two wells. In terms of the nuclear spin hamiltonian, 
this means that both the absolute value of the EFG and 
the quantization axis depend on the two-state variable. 
These quantities are denoted <^q and Uq in the ground 
state, and 0" and Ui in the excited level. 

Now we discuss how nuclear spins affect the polariza- 
tion echo that arises from a coherent superposition of 
the two tunnel states. In general, the quadrupolar part 
of the nuclear spin hamiltonian does not commute with 
H. Thus V leads to a dispersion of the two-level splitting 
when switching from Ug to Ui during the two pulses, and 
thus to a dephasing of the the echo signal. 

The density operator of a TLS involves four indepen- 
dent operators, e.g., the three Pauli matrices ax, (Jy, (Jz, 
and unity. Accordingly, the propagator is represented by 
a four-dimensional matrix p3| . Taking into account a 
nuclear spin / renders the dynamics significantly more 
complex, since we have to deal with D = 2(2/4- 1) quan- 
tum states corresponding to a density matrix with 
entries. It can be shown that a nuclear spin results in an 
overall factor of the echo amplitude 0, 



P(i,tw) = Po(t,iw)/(t,iw), (5) 
where /(t, t^) is determined by the nuclear spin energies 



/(t,iw) 



^e-'"7^(^7r2)e-»^*"/'i7^(^7rl) . (6) 



The bar indicates the ensemble average; time evolution 
for zero driving field is written in terms of the Liouville 
operator = (l/?i)[V, and the rotations TZ account 
for the two electric-field pulses of duration and Rabi 
frequency fi. Both C and TZ are superoperators that act 
on nuclear spin variables. The argument of TZ{6) is the 
"pulse area" 0. 

The first external-field pulse creates a coherent super- 
position of the two tunnel states, whereas the second 
pulse exchanges their phases. (In the Bloch spin pic- 
ture, this is related to rotations of the "polarization vec- 
tor" (a^).) The full density matrix is expressed through 
standard basis operators \vLia) (uj/?!, with z, j = 0, 1 and 
a, /5 = — /, ...,/. In the composite space of TLS and nu- 
clear spin variables, C and TZ are represented by tetrads 
Lpqrs and Rpqrs of dimension D^. Spelling out the ma- 
trices TZ and and taking the trace, we obtain the cor- 
rection factor 

f{t,tw)= ^ fafSjS cos[eoa - £if3)t/h - {eo^ - eis)t^/K\, 

(7) 

where nuclear spin energy levels (i.e. the eigenvalues of 
V) are denoted by Soa and eia and fap^yS depends on 
the matrix elements of TZ and CTz . (Details will be given 
elsewhere [Q.) Here we resort to a simple approxima- 
tion that is justified in various situations, such as short 
pulses or almost parallel quadrupolar quantization axes, 
and that is expected to grasp the essential physics in any 
case. Since the polarization echo occurs on a time scale 
much shorter than the waiting time, we may put t = t^. 
For a TLS initially in the ground state we the phase fac- 
tor, 



/(tw) 
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+ 1 4^, 



|Xa/3p|Xa5|^COs(w/35tw), (8) 



a, (3,5 



that depends on the overlaps 

Xa/3 = (uoa|ui/3) 
and the quadrupole spectrum, 



(9) 



(10) 



For the case where the system is initially in the excited 
state, we find a similar expression with £97 instead of 
ei^. Thus /(^w) describes the reduction of the whole 
echo signal. 

In the remainder of this paper, we discuss the reduc- 
tion factor (^) First we consider the situation where the 
quadrupole coupling is ineffective, such as a zero EFG, 
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parallel quantization axes Uq = Ui , or a very strong mag- 
netic field. Then the nuclear spin states corresponding 
to the tunnel levels are identical, juga) — |uia), for 
a = -/, ...,/, and Xa/3 = <5q/3, resulting in f{t^) = 1. 

In general, however, the EFG is finite and the 
quadrupolar hamiltonian is not the same for the two lev- 
els, resulting in a non-diagonal overlap matrix, Xq/3- The 
quadrupolar energy scale reads as 



Hluc 



4/(2/ - 1) 



(11) 



Typical values for the quadrupolar energy (('"Q corre- 
spond to frequencies of the order of tens of MHz and 
thus satisfy, for waiting times iw M sec, the inequality 
(^qiw ^ 27r that simplifies significantly the analysis. As a 
consequence, all terms involving different quadrupole lev- 
els 7 7^ ±/3 in (^) vanish. Yet note that the quadrupolar 
spectrum exhibits a degeneracy with respect to /? — > -~P 
which, in turn, is lifted by a magnetic field. 

In the limit of zero waiting time the cosines in (||) are 
equal to unity, and we have /(iw ^ 0) = 1. In the op- 
posite case of very long times and all degeneracies lifted, 
the terms with finite frequency Lu-y^ vanish, resulting in 



1 



oo 
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a,l3 



Xa/3 



1 



(12) 



We are interested in the intermediate regime of experi- 
mentally relevant waiting times that are of the order of 
IJ, sec. Thus we have to look for frequencies in the MHz 
range that satisfy the condition w^^jiw ~ i"- 

Both the overlaps Xa/i and the frequencies depend 
in an intricate manner on the relevant orientation of the 
three vectors e^, Uq, Ui and on the ratio of the Zeeman 
splitting and the quadrupolar energy. Here we discuss a 
few limiting cases where (^) simplifies significantly. The 
argument is developed for half-integer spin / = |i |: ■•• 
but easily generalized to integer /. 

First we consider the case of a weak magnetic field 
where 



(13) 



is small as compared to Uloq. Then the nuclear Zeeman 
splitting hujz lifts the degeneracy of the doublets ±/3 of 
the quadrupolar energy {(3 — i, /). Discarding rapidly 
oscillating terms ^ cos {ujqt^), and separating the weight 
factor 



and the time-dependent term, we obtain 

/(iw) = 1 - a + ^ b(}C0s{LJp-i3t^). 



(14) 



For small B we may neglect the magnetic-field depen- 
dence of the weight factors and treat the Zeeman term 
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FIG. 2: Magnetic-field dependence of the echo amplitude of 
BK7 jisj. The solid lines are calculated as explained in the 
text. 

as a perturbation. Starting from the eigenbasis of the 
quadrupolar energy, we thus diagonalizc the Zeeman 
term of V in each degenerate subspace ±/3. The splitting 
of each doublet ±/3 depends on the relative orientation 
of the quantization axes and ui through the cosine 
a; = • Ui. For /3 = |, /, the splitting is given by the 
projection of the magnetic field on the quadrupolar axis. 



cj/5,_;3 = x2[iuz {13 > 1/2), 
whereas for /? = i it reads as 



(15) 



- V^a;2 + (/+l/2)^(l-a;2)^z. (16) 
The average in (^^ is given by 



(...) 



dxp{x) (. 



(17) 



Though possible in principle, calculation of the normal- 
ized distribution p{x) is beyond the scope of the present 
paper. Eq. ( p^ shows oscillatory behavior with period 
^ IbJztvj, independent of the precise form of p{x). For 
the fit of the experimental data in Fig. 2, the best results 
are obtained with p(x) — 3x^, and bi = 2b3. Yet note 

^ ^ ' 2 2 

that the position of the first minimum of f(t^) occurs 
always at cjztw ~ O.Gtt and hardly depends on p{x). 

Now we turn to strong magnetic fields, luz 3> cjq, 
where the quadrupole energy may be treated as a pertur- 
bation with respect to the Zeeman splitting. When cal- 
culating the overlap matrix Xaf3 to lowest order in ujq/ujz 
and observing the normalization condition I Xa/3 P — 
1, we obtain 



1- /(<„) = 2 ^ Ixcpl 



(18) 
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Thus we find a variation /(^w) = 1 — const, x at high 
magnetic fields, as shown in Fig. 2. 

These theoretical findings (p^-p^ agree rather well 
with available data. We briefly discuss the most salient 
features. 

(i) Both the oscillatory behavior of the echo amplitude 
with small B and the saturation at higher fields have been 
observed for several multicomponent glasses and mixed 
crystals [|[ |^. In Fig. 2 we plot the echo amplitude mea- 
sured for the borosilicate glass BK7 as a function of the 
magnetic field. At small B the data show a few oscil- 
lations; at higher fields the amplitude increases strongly 
and would seem to saturate. The solid line for B < 100 
mT is calculated from Eq. (|lj) that, with t^j = 2/xsec, 
g = 1.8 and / = |, describes the first minimum unam- 
biguously. The increase at i? > 100 mT has been fitted 
with ®. 

(ii) The particular waiting time dependence of the echo 
amplitude through the product Btvj has been verified ex- 
perimentally in great detail, especially for KBr:CN 0. In 
Table I we give the nuclear spin parameters / and g and 
the measured and calculated positions of the first mini- 
mum in terms of the quantity iJmin^w Theoretical values 
are given for p{x) = 3a;^ as discussed above (p^). (The 
calculated values depend weakly on p{x) and bfj.) 

(iii) All systems showing the magnetic-field depen- 
dence contain nuclei with / > 1 and finite quadrupole 
moments (B in BK7; Br and N in KBr:CN; Al in Ba-Al- 
silicate). On the other hand, the only glass that shows 
no magnetic-field dependence (amorphous silicon oxide) 
[ p5[ , does not contain nuclear quadrupoles, since I = \ 
for 29Si and / = for and ^^Si. 



TABLE I: Nuclear spin parameters g and I. Experimental 
values for the first minimum of the echo amplitude. Calcu- 



lated values as 


explained 


in the text. 














(meas.) 


(calc.) 






/ 


g 


^miniw 


(10-^ Ts) 


borosilicate 


"B 


3/2 


1.8 


22 


21 


Al-Ba-silicate 


2^A1 


5/2 


1.44 


28 


22 


KBr:CN 


79,8lB^ 


3/2 


1.5 


30 


28 



These findings provide very strong evidence that the 
observed magnetic-field dependence arises from the dy- 
namic phase of the nuclear Zeeman splitting of quadrupo- 
lar levels. (Preliminary experiments on other systems 
would seem to confirm this statement The simpli- 

fications of the present theory may be at the origin of the 
discrepancies in Fig. 2. For example, real tunnel systems 
certainly involve more than one nuclear spin. This is ob- 
vious for glasses. In mixed crystals, the tunneling atom 



drags its dressing cloud; that's why bromine appears in 
Table I. 

It seems likely at this point that the quadrupole split- 
ting of tunneling levels not only influences the ampli- 
tude of polarization echoes, but has also consequences 
for other properties of glasses at very low temperatures. 
In particular, we expect that the magnetic field depen- 
dence of the dielectric susceptibility observed in several 
glasses is also caused by nuclear spins. 

In summary, we have proposed an explanation for the 
recently observed magnetic-field dependence of polariza- 
tion echoes in terms of a novel dephasing mechanism 
involving nuclear spins and quadrupole moments. The 
EFG (j>" and the corresponding axes u are not the same 
in the two minima of the atomic double-well potential; 
the quadrupolar energies gives rise to phase dispersion 
that reduces the echo signal. The oscillations at small 
fields result from the interference of the dynamical quan- 
tum phases of almost degenerate quadrupole levels. The 
strong increase of the echo amplitude at larger B is 
due to the alignment of the nuclear magnetic moments 
with respect to the magnetic field. For strong fields, 
pL^B 3> hujq^, we expect saturation at the value f ~ I. 
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